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"TO BLOCK OR CQVARY A CONCOMITANT VARIABLE: 
WHICH IS MORE POWERFUL?" 

The most widely used procedures to harness the power of 
a concomitant variable are block designs and ANCOVA. The 
decisions on whether to block or covary and how many blocks to be 
used if blocking is selected are often based on rules of thumb 
with little empirical support. The purpose of this study is to 
provide a scientific foundation on which to base such decisions. 

Monte Carlo generated data were analyzed using one-way 
ANOVA; two-block, four-block, and eight-block block designs; and 
ANCOVA. Resulting empirical powers were entered into a repeated 
measures four-way factorial design with three factors 
representing different experimental conditions and one factor 
representing the five procedures being compared. 

The results indicated that the correlation coefficient 
between the concomitant and dependent variables was the critical 
factor to influence the choice. One-way ANOVA was the best 
choice when there was no relationship, blocking was preferred 
when the correlation was low, and ANCOVA achieved the highest 
power when the correlation was high. Block designs and ANCOVA 
became more powerful and the optimal number of blocks increased 
as the correlation coefficient, the number of treatments, and the 
number of subjects per treatment increased. 
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TO BLOCK OR COVARY A CONCOMITANT VARIABLE: 
WHICH IS MORE POWERFUL? 

INTRODUCTION 

Educational experiments often involve assigning students to treatments. Traditional one-way analysis 
of variance can be used to analyze the differences among treatments. However, differences among 
students, such as, gender, socio-economic status, and level of ability, often mask or obscure the effects of 
a treatment (Kennedy & Bush, 1985; Kirk/ 1982). Nuisance variation due to such differences can be 
extracted from the error variance. By controlling the concomitant (nuisance) variable, researchers often 
reduce the background noise, increase the precision, and enhance the statistical power of a design (Bonett, 
.1982: Keppel, 1991; Maxwell & Delaney, 1984). The most widely used procedures to harness the power of 
a concomitant variable are block designs and the analysis of covariance. 

Statement of the Problem 

Evidence suggests that approximately 70% of published research in the behavioral sciences uses 
analysis of variance techniques (Glass & Hopkins, 1984). Many of these studies are analyses of either block 
or covariance designs. Whether to block or covary and how many blocks to be used if a block design is 
chosen become important questions and the answers may differ according to the experimental conditions. 
The purpose of this study is to determine which procedure should be used given a set of experimental 
conditions. 

Significance of the Study 

The decisions on whether to block or covary and how many blocks to be used if a block design is 
selected are often based on rules of thumb with little empirical support. An empirical study that can offer 
the scientific foundation on which to base such decisions is desirable. The results of this study should 
provide guidelines to help researchers decide the appropriate procedures to be used. 
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REVIEW OF THE RELATED LITERATURE 

Historical Review of the Problem 

By employing a concomitant variable, researchers <xn control the nuisance variance, reduce the error, 
increase the precision, and enhance the efficiency of an experimental design (Cochran & Cox, 1950; 
Cochran & Cox, 1957; Federer & Schlottfeldt, 1954; Fisher, 1937; Fisher, 1973a; Fisher, 1973b; Kennedy 
ana Bush 1985; Keppei, 1973; Keppel. 1991; Kirk, 1982; Undquist, 1953; Maxwell & Deianey, 1990; Myers, 
1979). In the two classic books. The Design of Experiments and Statistical Methods for Research Workers. 
Fisher (1937; 1973a) developed the analysis of variance of block designs and the analysis of covariance. 
He demonstrated that the precision of an experimental design could be improved by controlling a 
concomitant variable using the two analysis techniques. 

Undquist (1953) used the term, treatments-by-ievels design, which consists of more than one subject 
in a ceil, to differentiate it from the randomized complete block design, which consists of only one subject 
in a cell. The treatm nts-by-ievels design is also called the treatments-by-blocks design (Kennedy & Bush 
1985). Undquist recommended that the treatments-by-blocks design be used in preference to the analysis 
of covariance because (1) the treatments-by-blocks design required much less restrictive assumptions than 
the analysis of covariance, (2) the computational procedures were considerably simpler with the treatments- 
by-blocks design, and (3) the use of treatments-by-biocks design permitted a study on the simple effects 
of the treatments at any given block. 

Gouriay (1953) compared the analysis of covariance with the randomized complete block design in 
which blocks were formed by matching subjects on the concomitant variable. He recommended that the 
analysis of covariance be used in preference to the matching block technique; this view was shared by 
Greenberg (1953) in a similar study. 

Federer (1955) favored the block design over the analysis of covariance. He offered the following rule 
of thumb: "if the experimental variation cannot be controlled by stratification (blocking), then measure related 
variates and use covariance'' (p. 483-484). However, he also pointed out that "It may be more advantageous 
to use covariance than to use stratification, since fewer degrees of freedom are usually required to control 
the variation" (p. 484). 
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Cox (1957) developed the Apparent Imprecision measure and used it to compare the analysis of 
covariance with the randomized complete block design in which blocks were formed by ranking subjects 
on the concomitant variable. He found that the randomized complete block design was somewhat better 
than the analysis of covariance if the correlation coefficient was less than .6 while the analysis of covariance 
became appreciably better than the randomized complete block design when the correlation coefficient was 
.8 or more. He suggested that the analysis of covariance was preferable to the randomized complete block 
design only if the correlation coefficient between the concomitant and the dependent variable was at least 
.6. 

The most rigorous research on this topic was conducted by Feldt (1953). He used Cox's Apparent 
Imprecision measure to compare three experimental designs. The three experimental designs being 
compared were (1) stratification (blocking), (2) the analysis of covariance, and (3) the analysis of variance 
of difference scores. Feldt found the analysis of variance of difference scores was the least precise 
procedure; for p (correlation] < .4 block designs results in approximately equal or greater precision than 
the analysis of covariance; for p ^ .6 the advantage is in favor of the analysis of covariance; and for p < 
.2 and small values of the number of subjects neither the analysis of covariance nor block designs yields 
appreciably greater precision than the one-way analysis of variance. Feldt also provided a table for the 
optimal number of blocks to be used if block designs were selected. He concluded that the optimal number 
of blocks tended to be larger for (1 ) larger values of correlation coefficients, (2) larger numbers of subjects, 
and (3) smaller numbers of treatments. This study should be considered the classic study comparing block 
designs with ihe analysis of covariance; its findings have been most often quoted by textbooks in the area 
of experimental designs (e.g., Cook & Campbell, 1979; Dayton, 1970; Kennedy and Bush, 1985; Keppel, 
1991; Kirk, 1982; Myers, 1979). However, Feldt's concept of optimal blocking was not supported by 
Chuang's (1978) study, which found that block designs would become more powerful as the number of 
blocks increased and *the power of BLOCKING even at its maximum was slightly smaller than that of 
COVARIANCE" (p. 37). 




In a block design, subjects are usually grouped into blocks before the experiment according to the value 
of the concomitant variable. However, there are times that the value of the concomitant variable is not 
available before the experiment. When blocks are formed after the experiment, the block design is defined 
as a post-hoc block design. Keppel (1973) gave the following advantages of post-noc block designs over 
the analysis of covariance: (1) reduction in computational effort, (2) free from the stricter assumptions of the 
analysis of covariance, and (3) possibility of testing the treatments-by-blocks interaction. However, he also 
pointed out two disadvantages of post-hoc blocking: (1) the inability to calculate the within-groups mean 
square when cells have fewer than 2 subjects, and (2) the inability to adjust the treatment means for 
differences on the concomitant variable. 

Post-hoc blocking is popular because the value of the concomitant variable can be unknown before the 
experiment. Nevertheless, Myers (1979) pointed out the danger of abusing post-hoc block designs by 
demonstrating that reordering scores within each treatment does not change the treatment means but 
generally reduces the error variance, resulting in significant Fs which "merely reflect the reduction in error 
variance due to blocking rather than any variability due to treatments 11 (p. 155). However, he did not 
consider the loss of degrees of freedom with the block design. 

Bonett (1 982) compared post-hoc block designs with the analysis of covariance and offered the following 
rule: 

// the assumptions for each method can be satisfied and if the probability of a Type II error 
is of concern, the analysis of covariance will be preferred when the form of the regression 
equction is known but the magnitude of the correlation is known. Post-hoc blocking, on 
the other hand, will be preferred when the magnitude of the correlation is known, (p. 38) 

A study employing the Monte Carlo method and using statistical power as the criterion variable to 

compare block designs and the analysis of covariance was performed by Maxwell and Delaney (1 984). Their 

study was limited to two treatments. The procedures they compared were based on the following two 

dimensions: (1 ) the method of assignment and (2) the method of data analysis. Each of the two dimensions 

had three levels: the concomitant variable was (1) ignored, (2) categorized, and (3) continuous. This 

resulted in nine procedures being compared. Maxwell and Delaney (1984) favored the analysis of 

covariance over block designs. They argued that 
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the recommendation of most experimental design texts to consider the correlation between 
the dependent and concomitant variables in choosing the best technique for utilizing a 
concomitant variable is incorrect. Instead, the two factors that should be considered are 
whether scores on the concomitant variable are available for all subjects prior to assigning 
any subjects to treatment conditions and whether the relationship of the dependent and 
concomitant variables is linear, (p. 136) 

They also illustrated that the Apparent Imprecision measure, which was used in the Cox (1957) and Feidt 

(1958) studies, might provide a different perspective from statistical power, but, the Apparent Imprecision 

measure and statistical power are not independent. 

Summar/ 

While some research favored block designs, others preferred the analysis of covariance. Based on the 
historical review of the problem, it is summarized that "the relative merits of blocking and ANCOVA are more 
complicated, because neither is uniformly superior to the other - (Maxwell & Delaney, 1984, p. 136). It is 
likely that different procedures may be preferable to the others depending on the sets of experimental 
conditions. One consequence of applying block designs and the analysis of covariance, which has been 
often neglected in early research but frequently stressed in recent research, is the decrease of the probability 
of the Type 1! error, i.e., the increase in statistical power. 

Based on the review of the related literature, it is suggested that future research should examine the 
problem based on three dimensions: (1) how subjects are assigned, (2) how data are analyzed, ani (3) the 
distributions of and the relationship between the concomitant and the dependent variables (i.e., considering 
the assumptions of block designs and the analysis of covariance). The experimental conditions should 
include three factors: (1) the number of treatments, (2) the number of subjects per treatment, and (3) the 
magnitude of the relationship between the concomitant and the dependent variables. The criterion variables 
on which to base the comparison should be the statistical power, the Type I error (a), and the Apparent 
Imprecision measure. 

Justification of the Study 

This section provides the rationale for selecting statistical power as the criterion variable and using 
computer generated data 'o simulate the experiment. 
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Statistical Power as the Criterion Variable 

The expressions; "reduce error - , 'increase precision', 'enhance efficiency*, and 'maximize statistical 
power;' have frequently been used interchangeably to describe the objective of employing a concomitant 
variable in block designs and the anaJysis of variance (e.g., Bonett, 1982; Kennedy & Bush, 1985; Maxwell 
& Delaney, 1984). Among these expressions, the term "statistical power - is the most unambiguously 
understood one and is operationally defined in most statistical texts. 

The neglect of statistical power in research, textbooks, and curricula has been brought to the attention 
of the research community. As Cohen (1962; 1977; 1988; 1992) has stressed, one of the most pervasive 
threats to the validity of the statistical conclusions reached by behavioral research is low statistical power. 
T!^e investigation of statistical power in experiment designs has gained more and more popularity (Chase 
& Tucker, 1976; Sedlmeier & Gigerenzer, 1989). 
Computer Simulation 

This is an empirical study using the Monte Carlo method to simulate the experiment The Monte Carlo 
method has been used effectively in examining many properties of statistics (Harwell, Rubinstein, Hayes, 
& Olds, 1992; Shapiro, Wilk, & Chen, 1968; Wilcox, Charlin, & Thompson, 1986). Computer simulations have 
many advantages. "We can often simulate situations more readily on the computer than perform the 
corresponding experiments in real life"; "one can also easily vary parameters in computer experiments"; and 
"furthermore, the simulations tend to be very flexible in that a whole multitude of differing models can be 
simulated with relative ease with essentially the same computer code" (Jain, 1992, p. 2). Therefore, using 
a high speed computer to calculate statistical power based on empirical sampling is the most direct and 
effective way to answer the research questions of this study. 

PROCEDURES 

This study compared five analysis procedures under 48 sets of experimental conditions using empirical 
power as the criterion (dependent) variable. The five analysis procedures were the one-way analysis of 
variance; two-block, four-block, and eight-block block designs; and the analysis of covariar.ca. The 48 sets 
of experimental conditions were the combinations of four levels of the number of treatments (T; 2, 3, 4, 5), 
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three levels of the number of subjects per treatment (N; 8, 40, 72), and four levels of the correlation 
coefficient between the concomitant and the dependent variable (C; .00. .28, .56, .84). For each 
experimental condition, 1 ,000 sets of data were generated by the computer. Each set of data was analyzed 
by all five analysis procedures with the significance level (a) set at .05. The proportion of significant 
analyses was the empirical power; for example, if 600 out of the 1,000 analyses were significant, the 
empirical power would be .6. Each resulting empirical power was entered as one observation in each cell. 
The procedure to generate and analyze the 1,000 sets of data was repeated two more times for each of the 
48 experimental conditions. This resulted in a repeated measure four-way factorial design with three 
observations (i.e.. three empirical powers) In each of the 240 (5 X 48) cells. 



Statistical power is determined by three major factors: (1) the significance level, (2) the sample size, and 
(3) the effect size (Dayton, Schafer, & Rogers 1973; Hinkle, Wiersma, & Jurs, 1988; Upsey, 1990; Sawyer 
& Ball, 1981). Statistical power increases as the significance level, the sample size, and the effect size 
increases. In order to make the comparison of the five analysis procedures more meaningful, the one-way 
analysis of variance was treated as the control group by setting its power at .50. This also would allow the 
powers of the other procedures to increase or decrease as a function of the conditions. Therefore, the effect 
sizes which would achieve a .50 power for the one-way analysis of variance under given experimental 
conditions needed to be calculated before the experiment. Calculation of effect sizes was based on tables 
and formulae in Cohen's (1988) book. The effect size index is defined as 



where m, is the population mean of the ith treatment, m is the grand population mean, and k is the number 
of treatments. The effect size index was calculated by the following formula: 



Calculation of the Effect Size 



a 



where a is the population standard deviation and 




eric 



J u 



£ = 



n 



05 



\ 400 (n-1) ' 



where n is the number of subjects per treatment and n M is the required number of subjects per treatment 
to achieve a desired power when f is equal to .05. The value of n M can be obtained from Cohen's tables 
(p. 381 -389). For example, for three treatments with a = .05, the table shows that n M is 662. If the number 
of subjects per treatment is 8, then 



£ = 



\ 



662 



400 (8-1) 



=0 . 4862392 



When the means are equally spaced, the distance between the largest and the smallest mean can be 
calculated by the following formula: 



d = 2£ 



3 (jc-1) 
k+1 



In this case, 



d = 2 (0 .4862392), 



3 (3-1) 
N 3+1 



=1. 1910379 



in order to equally space the means, this number should be divided by two: 



1 . 1910379 



=0.5955189 



Thus, to achieve a .5 statistical power for one-way analysis of variance when a =.05, n = 8, and k = 3, the 
population mean of the first treatment should be set at m„ the second at m, + 0.5955, and the third at m, 
+ 1.1910. Since the results of the analyses would be the same for any value of m,, the population mean 
of the first treatment was always set at 0 in this experiment The following is the table of the calculated 
effect sizes to achieve a statistical power of .50 for the one-way analyses of variance under given 
experimental conditions. 
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Table 1 



Calculated affect Size 





Effect sizes for N=8 


effect sizes for M=40 


effect sizes for Ms 72 




1 .0481 


0.U40 


0,3291 | 


T=3 


0.5955, 1.1910 


0.2523, 0.5046 


0.1870, 0.3740 


T=4 


0.4060, 0.8121, 1.2181 


0.1720, 0.3440, 0.5161 


0.1275, 0.2550, 0.3825 


T=5 


0.3030^ 0.6059, 0.9089, 1 .2118 


9.1284^0,2567, 0.3851^0.5134 


0.0951^ 0.1903, 0.2354. 0.3805 J 



feneration and Analyses of the Data 
The generation and the analyses of the data were accomplished using a computer simulation system 
running on an IBM 3090/400E mainframe computer. Bivariate correlated data were generated using the SAS 
commands provided by Clark and Woodward (1992). These commands generated random data from a 
bivariate (the concomitant and the dependent variables) normal distribution with means of 0, variances of 
1, and a user-specified correlation coefficient. Random samples were generated separately for each 
treatment. Only the means of the dependent variable were transformed based on the calculated effect sizes 
in Tabie 1, while the other parameters remained unchanged. Data in each treatment were grouped into 2, 
4, and 8 blocks by their ranks on the concomitant variable. For example, to group 40 subjects into 4 blocks, 
the top 10 ranked subjects would be in the first block, the 11-20 ranked subjects would be in the second 
block, the 21-30 ranked subjects would be in the third block, and the 31-40 ranked subjects would be in the 
fourth block. 

The computer simulation system included one executable file and two SAS programs (International 
Business Machines, 1988a; Internationa] Business Machines, 1988b; SAS Institute Inc., 1990a; SAS Institute 
Inc., 1990b). For each of the 48 experimental conditions, the executable file ran the first SAS program 1 ,000 
times, then ran the second SAS program. The computer code for the condition of the number of treatments 
(T) equal to 5, the number of subjects per treatment (N) equal to 72, and the correlation coefficient (C) equal 
to .84 are listed in Appendix A, The first SAS program generated a set of data, analyzed that set of data 
with the five analysis procedures being compared, and output the results of the analyses to a data file. After 
the first SAS program ran 1,000 times, the data file contained 1,000 records of the results of the analyses. 

ERLC J * 



10 

The second SAS program calculated the empirical power based on the 1,000 records, in order to obtain 
three observations per ceil, the executable file ran three times for each of the 48 experimental conditions. 
Totally, there were 144,000 (1,000 X 3 X 46) sets of data generated and 720,000 (5 X 144,000) analyses 
conducted. 

A seed must be provided to generate random data using SAS random functions. The values of seeds 
can be any integer ranging from 1 to 2 31 - 2 (i.e., 2,147,483,646). Users can let the computer clock set an 
initial seed by specifying a value of 0. Using the computer clock to generate random data was tested and 
found to have three problems: (1) The computer clock may not increment enough to generate different data; 
(2) The computer clock may generate repeated or patterned data; and (3) the program may not be executed 
because the computer clock generates invalid seeds. Therefore, positive seeds were used instead of the 
computer clock. Using positive seeds also makes the experiment repilcable. in order to systematically and 
representatively employ the seeds, the minimum seed value, 1, was useu as the first initial seed; it was 
incremented by 2,147,483 for each run of the first. SAS program, and by 14,913 far each run of the 
executable file. The two incremented values were obtained by dividing the maximum seed value by 1,000 
and 144,000. Thus, all seeds were equally spaced in the range between the minimum and maximum seed 
values. The initial seeds used for each executable file are listed in Appendix B, 

Hypotheses 

This study tested the following null hypotheses using empirical power as the dependent variable. If a 
null hypothesis was rejected, its follow-ups were conducted, and the null hypotheses following it was 
ignored. 

HO r - There will be no significant differences at the .05 level for the four-way interaction. 
HO^ There will be no significant differences at the .05 level for the three-way interaction. 
HOy There will be no significant differences at the .05 level for the two-way interaction. 
H0 4 : There will be no significant differences at the .05 level for the main effect 

RESULTS 

The resulting power values are listed in Appendix C. The raw data were analyzed to test the null 
hypotheses. The results of the analyses are summarized in Table 2. 
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Table 2 
ANOVA Summary 



Tests of hypotheses using S(T*N*C) as the error term Tests of hypotheses using P*S(T*M*C) as the error term 



Source 


OF 


SS 


Mean Square 


F Value 


Pr > F 


Source 


OF 


SS 


Mean Square 


F Value 


Pr > F 


T 


3 


0.167672 


0.05589083 


116.31 


0.0001 


? 


4 


2.226638 


0.55667191 


12117.43 


0.0001 


N 


2 


0.050754 


0.02537705 


52.81 


0.0001 


T*P 


12 


0.073263 


0.00610521 


132.90 


0.0001 


C 


3 


8.853233 


2.95107773 


6141 .08 


0.0001 


N*P 


8 


0.026325 


0.00329068 


71.63 


0.0001 


T*N 


6 


0.010605 


0.00176756 


3.68 


0.0025 


C*P 


12 


2.820703 


0.23505859 


5116.69 


0.0001 


T*C 


9 


0.103833 


0.01153697 


24.01 


0.0001 


T*N*P 


24 


0.005493 


0.00022889 


4.98 


0.0001 


N*C 


6 


0.002530 


0.00042162 


0.88 


0.5146 


T*C*P 


36 


0.051654 


0.00143483 


31.23 


0.0001 


T*N*C 


18 


0.004716 


0.00026198 


0.55 


0.9285 


N*C*P 


24 


0.004126 


0.00017191 


3.74 


0.0001 


S<T*N*C) 


96 


0.046133 


0.00048055 






T*M*C*P 


72 


0.002742 


0.00003808 


0.33 


0.8336 














P*S(T*H*C) 


384 


0.017641 


0.00004594 


















Total 


719 


14.468110 









Note . T: Number cf Treatments, N: Number of Subjects per Treatment, C: Correlation Coefficient, P: Procedure, 
S: Oata Set. 



The four-way interaction (T*N*C*P) was not significant while the three-way interactions of T*N*P, T*C*P, 
and N*C*P were significant. The cell means for the significant three-way interactions are listed in Tables 
3, 4, and 5 respectively. Simple simple effects were tested at the significance level of .01. Non-significant 
simple simple effects are indicated by "NS" in the last cells of the corresponding rows and columns. Tukey's 
Honest Significant Difference (HSD) values are provided in the table notes for multiple comparisons. The 
cell means for the other combinations are listed in Appendix D. 

Table 3 



• Means for the Interaction of the Nunber of Treatments, the N umber of Subjects per Treatment, and the Procedure 
(T*H»Q) 





ANOVA 


TWO -BLOCK 


FOUR -BLOCK 


SIGHT -BLOCK 


ANCOVA 


T2 


N08 


.498 


.561 


.580 


.546 


.628 


N40 


.500 


.575 


.604 


.613 


.663 


N72 


.506 


.577 


.606 


.615 


.663 


T3 


M08 


.497 


.583 


.606 


.599 


.646 


N40 


.501 


.593 


.627 


.639 


.666 


N72 


.502 


.593 


.628 


.639 


.666 


T4 


nob 


.501 


.601 


.633 


.630 


.653 


N40 


.502 


.609 


.645 


.657 


.668 


N72 


.501 


.613 


.650 


.663 


.674 


T5 


N08 


.503 


.615 


.649 


.648 


.666 


N40 


.504 


.626 


.663 


.677 


.674 


N72 


.498 NS 


.620 


.656 


.668 


.671 



Note . P3T*N = .008, T*NS)P -.015, NS: Non-significant. 
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Table 4 

Means for the Interaction of The Mumper of Treatments, the Correlation Coeff*c ; ent. and the Procedure (T»C*P) 





ANOVA 


TWO-BLOCK 


FOUR-BLOCK 


SIGHT-SLOCK 


AMCOVA 


T2 


COO 


.504 


.502 


.498 


.431 


.490 


C28 


.502 


.518 


.523 


.511 


.523 


CS6' 


.496 


.576 


.605 


.601 


.645 


C84 


.503 


.689 


.759 


.774 


.947 


T3 


COO 


.505 


.505 


.501 


.495 


.496 


C23 


.497 


.519 


.525 


.522 


.523 


C56 


.502 


.601 


.636 


.643 


.653 


C84 


.496 


.735 


.821 


.844 


.960 


T4 


COO 


.503 


.503 


.504 | .496 


.498 MS 


C28 


.496 


.527 


.535 1 .534 


.526 


CS6 


.498 


.623 


.663 


.676 


.677 


C84 


.508 


.778 


.868 


.394 


.967 


T5 


COO 


.505 


.503 


.504 


.499 


.497 MS 


C28 


.501 


.534 


.544 


.545 


.526 


C56 


.507 


.640 


.686 


.699 


.683 


C84 


.494 MS 


.804 


.890 


.916 


.975 



Mote . HSO: P3T*C=.009, T*caP=.019, MS: Mon-signif icant. 
Table 5 

Means for the Interaction of the Mumber of Subjects per Treatment, the Correlation Coefficient, and the 
Procedure_(N»C*P) 





ANOVA 


TWO- BLOCK 


FOUR-BLOCK 


EIGHT-BLOCK 


AMCOVA 


M08 


COO 


.502 


.499 


.494 


.469 


.478 


C28 


.500 


.521 


.525 


.506 


.512 


C56 


.498 


.602 


.634 


.624 


.653 


C84 


.499 


.739 


.816 


.324 


.955 


M40 


COO 


.505 


.505 


.505 


.505 


.503 MS 


C28 


.502 


.527 


.537 


.541 


.534 


C56 


.502 


.614 


.652 


.667 


.670 


C84 


.498 


.759 


.844 


.374 


.964 


M72 


COO 


.505 


.506 


.505 


.504 


.505 MS 


C28 


.495 


.525 


.534 


.537 


.530 


C56 


.502 


.615 


.657 


.673 


.670 


C84 


.503 MS 


.758 


.844 


.872 


.968 



Note . HSO: P3N*C=.008 r M*C3P=.015, MS: Mon-signi f icant 
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Comparison of Analysis Procedures 
When the correlation coefficient between the concomitant and dependent variables was equal to .00, 
the one-way ANOVA was as powerful as or more powerful than blocking and covariance. The difference 
was significant when the number of treatments and the number of subjects per treatment were small. As 
the correlation coefficient increased, the power of employing a concomitant variable became more and more 
significant— from 50 to over 90 percent. With the correlation coefficient equal to .28, the optimal blocking 
procedure was as powerful as or slightly more powerful than covariance. Covariance was favored when the 
correlation was moderate or high. However, with moderate correlation, blocking could be as powerful as 
or slightly more powerful than covariance when the number of treatments and the number of subjects per 
treatment were large. 

Comparison of Conditions 
The ranks of the power values of the 48 sets of conditions for all and each of the analysis procedures 
are listed in Appendix E. The pattern of the ranks showed that the three blocking procedures and the 
ANCOVA became more powerful as (1) the correlation coefficient, (2) the number of treatments, and (3) the 
number of subjects per treatment increased. Among the three factors t the correlation coefficient was 
dominant The ranks for the one-way ANOVA showed a random pattern with all power values being 
approximately .50, which provided evidence that the power of the one-way ANOVA had been successfully 
controlled. 

Optimal Number of Blocks 
The results did not provide the specific optimal number of blocks to be used under each experimental 
condition because the experiment did not include ail possible numbers of blocks. The results did indicate 
that the optimal number of blocks increased as the correlation coefficient, the number of subjects per 
treatment, and the number of treatments increased. 

SUMMARY 

Different procedures should be used depending on the set of experimental conditions. The correlation 
coefficient between the concomitant and the dependent variable was the critical factor that should influence 
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the choice. The one-way ANOVA was the best choice when the correlation was zero, block designs were 
preferred when the correlation was low, and the analysis of covariance achieved the highest power when 
the correlation was high. With moderate correlation, block designs should be selected only when the 
number of treatments and the number of subjects per treatment were large; otherwise, the analysis of 
covariance should be used. Block designs and the analysis of covariance became more powerful and the 
optimal number of blocks for a block design increased as xhe correlation coefficient the number of 
treatments, and the number of subjects per treatment increased. 

Discussion 

The levels of the experimental conditions were chosen to be equally spaced and to be representative 
of real world situations. The four levels of the number of treatments represented the most commonly used 
numbers of treatments; the three levels of the number of subjects per treatment represented small, medium, 
and large sample sizes; and the four levels of the correlation coefficient represented zero, low, moderate, 
and high correlations. The results of the study provided a guide to help researchers decide the appropriate 
procedures to be used under different experimental conditions. 

This study had the following characteristics: 

1 . It controlled the power of the one-way ANOVA to prevent ceiling and floor effects; this also made the 
comparisons more meaningful as the one-way ANOVA was treated as the control group. 

2. Unlike most of the Monte Carlo studies that had only one observation in each cell and provided only 
descriptive statistics, this study had multiple observations in each cell and provided inferential in addition 
to descriptive results. 

3. Basing the sample size on the number of subjects per treatment in?«;ead of the total number of subjects 
made the interpretation of the results more meaningful. 

4. The computer simulation system consisted only of SAS programs, which were much shorter and more 
understandable than equivalent programs written in Fortran or other programming languages. 

One limitation of this study was that it did not include all possible blocking procedures. Since the results 
showed that the optimal number of blocks increased as the sample size increased, block designs could 
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become more powerful if other blocking procedures with more blocks were used. The primary disadvantage 
of using SAS rather than a programming language for Monte Carlo simulation is that SAS uses more 
computer CPU time. However, this disadvantage can be overcome by using a high speed computer such 
as the IBM 3090/400E computer used in this study. 

Recommendations for Future Research 
Several pilot studies which examined the parameters and distributions of the mean, variance, and 
correlation coefficient were conducted before the experiment, and the resulting sampling distributions of the 
statistics were checked after this experiment. The inspection found that the computer simulation system 
generated data that met predetermined specification, Furthermore, before the experiment the power of the 
one-way ANOVA was controlled at .50. The resulting empirical powers of the one-way ANOVA had a mean 
of .50 and a variance of .0001; aJso the mean squares for S(T*N*C) and P*S(T*N*C) were .00048055 and 
.00004594 respectively and the pooled mean square error was .00013286-aiI supporting the precision of the 
data generation procedures. Since the computer simulation system is able to generate accurate data and 
examine the problems effectively, it is recommended that future research adapt the system to examine 
related problems. 

This study does not include the treatments-by-biocks interaction in the block designs since the 
interaction does not exist in the population. Future studies can examine the effects of including the 
interaction using essentially the same computer codes t or, by varying the parameters of the population, 
examine the effects of including and excluding the interaction when the interaction exists in the population. 
The optimal number of blocks for a block design could be investigated by including other feasible blocking 
schemes such as 5-, 10-, 20-, and 40-biock block designs for the condition of 40 subjects per treatment 

The greatest contribution of this study may not be the specific results reported here, but the potential 
for examining many other situations. The computer simulation system developed for this study could be 
modified easily for a multitude of other studies. For example, it could be used to investigate other criteria 
such as Type I errors, examine other levels of the experimental conditions, or test other blocking methods 
in addition to the post-hoc blocking used in this study. 
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Appendix A 
COMPUTER CODES 



Executable File 

/* */ 
ADDRESS COMMAND 
"ERASE PVALUE DATA A " 
NUMERIC DIGITS 10 
TIME = 1 

DO WHILE TIME < 1001 

SEED a 2132560 + (TIME -1) * 2147483 

"EXEC 10 1 DISKW" NEWSEED DATA A 

" ( STRING " SEED 

"EXEC SAS T57284" 

"ERASE NEWSEED DATA A" 

TIME=*TIME+1 

END 

"EXEC SAS T57284P" 

First SAS Program (TS7284 SAS A) 

CMS FILEDEF INDATA DISK NEWSEED DATA 

A; 

CMS FILEDEF PVALUE DISK PVALUE DATA A 
(LRECL 210 BLKSIZE 210 RECFM FBS ; 
CMS FILEDEF SASLIST DISK T57284 
LISTING A; 

DATA BIVNORM (DROP~I); 
INFILE INDATA; 

INPUT SEED; 
DO I»l TO 72; 

GROUP- 1; 

X^RANNOR ( SEED ) ; 

Y* . 84*X+SQRT ( 1- . 84**2 ) *RANNOR( SEED) ; 
OUTPUT; 
END; 

DO I»l TO 72; 
GROUP-2; 
X=RANNOR ( SEED ) ; 

Y=».34*X+SQRT( 1-. 84**2 )*RANNOR( SEED) ; 
Y=0.Q951+Y; 



OUTPUT ; 
END ; 

DO 1=1 TO 72; 
GROUP = 3; 
X=RANNOR ( SEED ) ; 

Y= . 84*X+SQRT ( 1- . 84**2 ) *RANNOR ( SEED ) ; 
Y=0. 1903+Y; 
OUTPUT ; 
END; 

DO I»l TO 72; 
GROUP-4 ; 
X-RANNOR { SEED ) ; 

Y= • S4*X+SQRT ( 1- . 84**2 ) *RANNOR ( SEED ) ; 
Y=0.2854+Y; 
OUTPUT ; 
END; 

DO 1*1 TO 72; 
GROUP* 5 ; 
X 3 RANNOR ( SEED ) ; 

Y-.84*X+SQRT( 1-. 84**2) *RANNOR ( SEED ) ; 
Y~0.3805-t-Y; 
OUTPUT; 
END; 
PROC SORT; 

BY GROUP X; 
DATA BIVNORM; 
SET BIVNORM; 

BN»MOD<_N_,72) ; IF BN=*0 THEN 
BN=72; 

IF BN<=*36 THEN B2»l; ELSE B2»2; 

IF BN<=18 THEN B4-1; ELSE IF 

BN<=36 THEN B4-2; 

ELSE IF BN<- 54 THEN B4*3; ELSE 

B4»4; 

IF BN<=9 THEN B8*l; ELSE IF BN<=»18 
THEN B8=2; 
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ELSE I? 3N<=2 7 THEN 38=3; ELSE IF 
3N<=36 THEN 38=4; 

ELSE I? BN<=4 5 THEN B8=5 ; ELSE IF 
3N<=54 THEN B3=6; 

ELSE IF BN< = 63 THEN 38=* 7 ; ELSE 
B8=8; 

PROC PRINT; 

PROC CORK D AT A= B I VNORM ; 

VAR X Y ; 

BY GROUF; 
PROC GLM; 

CLASS GROUP ; 

MODEL Y=GROU?/SS3; 
PROC GLM; 

CLASS GROUP 32; 

MODEL Y=GROU? B2/SS3; 
PROC GLM; , 

CLASS GROUP 34; 

MODEL Y=GROUP 34/SS3; 
PROC GLM; 

CLASS GROUP 38; 

MODEL Y=GROU? 38/SS3; 
PROC GLM; 

CLASS GROUP ; 

MODEL Y=GROUP X/SS3; 
DATA; 

INFILE SASLIST; 

INPUT WORD] $ WORD 2 $ @; 

FILE P VALUE MOD; 

IF WORD1 » 'X' AND WORD 2 ='72' 

THEN DO; 

INFUT MEAN STDDEV; 

PUT MEAN 6.4 STDDEV 6.4 @; 

INPUT Y $ N MEAN STDDEV; 

PUT MEAN 6.4 STDDEV 6.4 @; 

END; 

ELSE IF WORDl^X" AND WORD 2 =» 
'1.00000' THEN DO; 
INPUT CORK; 
PUT CORK 6.4 @; 
END ; 

ELSE IF WORD 1='* GROUP M AND WORD 2 » 
'4' THEN DO; 
INPUT SS MS F PR; 
PUT PR 6.4 @; 

INPUT BLOCX $ DF SS MS F PR; 
PUT PR 6.4 6; 
END ; 

Second SAS Program (T57284P SAS A) 

CMS FILEDEF INDATA DISK PVALUE DATA 
A ; 

DATA PVALUE; 
INFILE INDATA; 

INPUT (G1XMEAN G1XSD G1YMEAN GlYSD 
G1CORR 

G2XMEAN G2XSD G2YMEAN G2YSD 
G2 CORR 

G3XMEAN G3XSD G3YMEAN G3YSD 
G3CORR 

G4XMEAN G4XSD G4YMEAN G4YSD 
G4CORR 



GSXMEAtf GSXSD G5YMEAN G5YSD 
G5CORR 

GROUP 13 BLOCK1B GROUP 2 3 
3LOCK2B GROUP4B BLOCK4B 
GROUP 8B 3LOCK8B GROUPANC 
3LOCKANC) (35* 6.4); 

TOTAL=0 ; 

G1BSG=0; 

B1BSG=0; 

G23SG=0; 

B2BSG=0; 

G4BSG=0? 

B4BSG=0; 

G8BSG=0; 

B8BSG=0; 

GANCSG^O ; 

BANCSG=0 ; 

TOTAL =1 ; 

IF GROUP 13 <= 0.05 THEN G13SG=1; 

IF BLOCK1B <= 0.05 THEN B1BSG=1; 

IF GROUP 2 B <= 0.05 THEN G2BSG=1; 

IF BLOCK23 <= 0.05 THEN 323SG=1; 

IF GROUP4B <= 0.05 THEN G4BSG=1; 

IF BLOCK4B <= 0.05 THEN B4BSG=1; • 

IF GROUP 8B <= 0.05 THEN G8BSG=1; 

IF BLOCK8B <= 0.05 THEN B8BSG=1; 

IF GROUPANC <= 0.05 THEN 
GANCSG=1; 

IF BLOCKANC <= 0.05 THEN 
BANCSG-1; 
PROC FREQ; 

TABLE G1BSG « BANCSG; 
PROC SUMMARY D AT A= PVALUE ; 

VAR G1XMEAN — BANCSG; 
. OUTPUT OUT =« DESCRIPT; 
PROC PRINT DATA=DESCRIPT; 
PROC UNIVARIATE DATA=P VALUE PLOT 
NORMAL; 

VAR G1XMEAN — BLOCKANC; 
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Appendix B 
INITIAL SEEDS 



SEED 



CONDITION 



SEED 



CONDITION* 



SEED 



CONDITION* 



1 

14914 

29827 

4/ "40 

59653 

74566 

89479 

104392 

119305 

134218 

149131 

164044 

178957 

193870 

208783 

223696 

238609 

253522 

26S435 

283348 

298261 

313174 

328087 

343000 

357913 

372826 

387739 

402652 

417565 

432478 

447391 

462304 

477217 

492130 

507043 

521956 

536869 

551782 

566695 

581608 

596521 

611434 

6263^7 

641260 

656173 

671086 

635999 

700912 



T2N08C00DS1 

T2N08C28DS4 

T2N08C56DS7 

T2N08C34DS10 

T2N40C00DS13 

f2N40C28DS16 

T2N40C56DS19 

T2N40C34DS22 

T2N72C00DS25 

T2N72C28DS28 

T2N72C56DS31 

T2N72C34DS34 

T3N08C00DS37 

T3N08C28DS40 

T3N08C56DS43 

T3N08C84DS46 

T3N40C0ODS49 

T3N40C28DSS2 

T3N40C56DS55 

T3N40C84DS58 

T3N72C00DS61 

T3N72C28DS64 

T3N72CS6DS67 

T3N72C84DS70 

T4N08C00DS73 

T4N08C28DS76 

T4N08C56DS79 

T4N08C84DS82 

T4N40C.00DS85 

: *N40C28DS88 

T4N40C56DS91 

T4N40C84DS94 

T4N72C00DS97 

T4N72C28DS100 

T4N72C56DS103 

T4N72C8 4DS106 

T5N08C00DS109 

T5N08C28DS112 

T5N08C56DS115 

T5N08C34DS118 

T5N40C00DS121 

T5N40C28DS124 

T5N40C56DS127 

T5N40C84DS130 

T5N72C00DS133 

T5N72C28DS136 

T5N72C56DS139 

T5N72C34DS142 



715825 

730738 

745651 

760564 

775477 

790390 

805303 

820216 

835129 

850042 

864955 

879863 

894781 

909694 

924607 

939520 

954433 

969346 

984259 

999172 

1014085 

1028998 

1043911 

1058824 

1073737 

1088650 

1103563 

1118476 

1133389 

1148302 

1163215 

1178128 

1193041 

1207954 

1222867 

1237780 

1252693 

1267606 

1282519 

1297432 

1312345 

1327258 

1342171 

1357084 

1371997 

1386910 

1401823 

1416736 



T2N08C00DS2 

T2N08C28DS5 

T2N08C56DS8 

T2N08C84DS11 

T2N40C00DS14 

T2N40C28DS17 

T2N40C56DS20 

T2N40C34DS23 

T2N72C0CDS26 

T2N72C28DS29 

T2N72C56DS32 

T2N72C84DS35 

T3N08C00DS38 

T3N08C28DS41 

T3N08C56DS44 

T3N08C84DS47 

T3N40C00DS50 

T3N40C28DS53 

T3N40C56DS56 

T3N40C84DS59 

T3N72C00DS62 

T3N72C28DS65 

T3N72C56DS68 

T3N72C34DS71 

T4N03C00DS74 

T4N08C28DS77 

T4N08C56DS80 

T4N08C34DS83 

T4N40C00DS86 

T4N40C28DS89 

T4N40C5 6DS92 

T4N40C34DS95 

T4N72C00DS98 

T4N72C28DS101 

T4N72C56DS104 

T4N72C84DS107 

T5N08C00DS110 

T5N08C28DS113 

T5N08C56DS116 

T5N08C84DS119 

T5N40C00DS122 

T5N40C28DS125 

T5N40C56DS128 

T5N40C84DS131 

T5N72C00DS134 

T5N72C28DS137 

T5N72CS6DS140 

T5N72C84DS143 



1431649 
1446562 
1461475 
1476388 
1491301 
1506214 
1521127 
1536040 
1550953 
1565866 
1530779 
1595692 
1610605 
1625518 
1640431 
1655344 
1670257 
1685170 
1700083 
1714596 
1729909 
1744822 
1759735 
1774648 
1789561 
1804474 
1819387 
1834300 
1849213 
1864126 
1879039 
1893952 
1908865 
1923778 
1938691 
1953604 
1968517 
1983430 
1998343 
2013256 
2028169 
2043082 
2057995 
2072908 
2087821 
2102734 
2117647 
2132560 



T2N08C00DS3 

T2N08C28DS6 

T2N08C56DS9 

T2N08C84DS12 

T2N40C00DS15 

T2N40C28DS18 

T2N40C56DS21 

T2N40C84DS24 

T2N72C00DS27 

T2N72C28DS30 

T2N72C56DS33 

T2N72C84DS36 

T3N08C00DS39 

T3N08C23DS42 

T3N08C56DS45 

T3N08C84DS48 

T3N40C00DS51 

T3N40C28DS54 

T3N40C56DS57 

T3N40C84DS60 

T3N72C00DS63 

T3N72C28DS66 

T3N72C56DS69 

T3N72C84DS72 

T4N08C00DS7 5 

T4N08C28DS7 8 

T4N08C56DS81 

T4N08C84DS84 

T4N40C00DS87 

T4N40C28DS90 

T4N40C56DS93 

T4N40C84DS96 

T4N72C00DS99 

T4N72C28DS102 

T4N72CS6DS10S 

T4N72C84DS108 

T5N08C00DS111 

T5N08C28DS114 

T5N08C56DS117 

T5N08C34DS120 

T5N40C00DS123 

T5N40C28DS126 

T5N40C56DS129 

T5N40C84DS132 

T5N72C00DS135 

T5N72C28DS138 

T5N72CS6DS141 

T5N72C84DS144 



* T: the number of treatments, N: the number of subjects per treatment, 
C: the correlation coefficient, DS: the data set. 
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C: the correlation coefficient, OS: the data set. 
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Mote. T: the numoer or treatments, N : the nunoer or suojects per treatment, 
C: the correlation coefficient, DS: the data set. 
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Mote. T: the number of treatments, M: tne numoer of suojects per treatment, 
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0.507 


T2M72C00 


0.510 


T4N08COO 


0.497 


72N72C00 


0.511 


75M40C00 


0.511 


72K72C00 


0.509 


74N72C00 


0.503 


75N40CQ0 


0.510 


T2M40C84 


0.498 


73M08C28 


0.514 


73N08C28 


0.516 


75N40C00 


0.513 


73H08C28 


0.512 


73N72C28 


0.518 


T3N0SC56 


0.498 


73M72C28 


0.516 


73N72C28 


0.525 


T4NC5C28 


0.522 


74W08C23 


0.518 


74N08C28 


0.521 


T4N72C56 


0.498 


72N40C28 


0.521 


72N40C28 


0.529 


73H72C28 


0.529 


75N08C28 


0.518 


72N40C28 


0.523 


T5N72C00 


0.498 


72M72C28 


0.522 


72N72C28 


0.529 


T2N72C28 


0.531 


74N72C28 


0.527 


T2N72C28 


0.524 


T5N72C28 


0.498 


74N40C28 


0.523 


74N08C28 


0.532 


72N40C23 


0.532 


75N72C28 


0.528 


T4N40C28 


0.525 


T2H40C56 


0.499 


T3M40C28 


0,526 


73N40C28 


G.533 


T5N08C28 


0.533 


73K72C28 


0.530 


74N72C28 


0.525 


T2S72C56 


0.499 


74N08C28 


0.529 


74M40C2B 


0.536 


73N40C28 


0.537 


74N40C28 


0.532 


T5M08C28 


0.525 


T3N40C56 


0.500 


74M72C28 


0.529 


74M72C28 


0.537 


74N72C28 


0.539 


75N40C28 


0.532 


73N40C28 


0.529 


T2M08C84 


0.501 


75JJ08C28 


0.532 


75M08C28 


0.540 


74N40C28 


0.540 


72N40C28 


0.533 


75N72C28 


0.530 


T2S40C28 


0.501 


75M72C28 


0.532 


75N72C28 


0.344 


75N72C28 


0.547 


72M72C28 


0.535 


75N40C28 


0.535 


T2N72C28 


0.501 


75M40C28 


0.537 


75N40C28 


0.549 


72M08C56 


0.554 


73N40C28 


0.540 


T2N08CS6 


0.564 


74M08C28 


0.501 


72JJ08C56 


0.565 


T2N08C56 


0.589 


75N40C28 


0.555 


T2H08C56 


0.623 


T3N08C56 


0.590 


74M40C56 


0.5Q2 


T2H72"' 


0.579 


72N40C56 


0.613 


73N08C56 


0.607 


73N08C56 


0.634 


72N40C56 


0.595 


T5N08C28 


0.502 


72K40Ljo 


0.583 


72N72C56 


0.614 


72N40C56 


0.622 


72N72C56 


0.655 


72H72CS6 


0.595 


72N40C0G 


0.503 


73M08C56 


0,591 


73N08C56 


0.617 


72K72C56 


0.627 


72N40C56 


0.657 


73N40C56 


0.613 


T3N40C00 


0.503 


73M40C56 


0.602 


73N40C56 


0.641 


73N40C56 


0.656 


T3N72C56 


0.653 


T4X08C56 


0.617 


73N72C30 


0.504 


73M72C56 


0.610 


73M72C56 


0.648 


74N08C56 


0.657 


74N08C56 


0.666 


T3M72C56 


0.618 


T3M72C84 


0.504 


74N08C56 


0.616 


i»*^08C56 


0.65 ; 


73N72C56 


0.665 


i 0N40Cjo 


rt LCI 

U.667 


74N40C56 


0.628 


75M40C28 


0.504 


T4N40C56 


0.625 


74W40C56 


0.662 


74N40C56 


0.679 


T4M40CS6 


0.672 


74N72C56 


0.637 


T5N72C56 


0.504 


74M72C56 


0.626 


74H72CS6 


0.677 


75N08C56 


0.680 


75N72C56 


0.675 


T5M08C56 


0.638 


T2M08C28 


0.505 


75M08CS6 


0.633 


T5M08C56 


0.678 


74N72C56 


0.693 


T5N40C56 


0.685 


T5M72C56 


0.643 


74W40C00 


0.505 


75M72C56 


G.642 


75M40C56 


0.690 


75N72C56 


C.706 


75N08C56 


0.689 


75N40C56 


0.648 


74M40C84 


0.505 


75M40C56 


0.645 


75M72C56 


0.690 


75N40C56 


0.710 


T4N72C56 


0.692 


72N08CS4 


0.717 


T4M72C0O 


0.506 


72M08C84 


0.680 


72N08C84 


0.740 


72N08C84 


0.730 


T2N08C84 


0.934 


72N40CS4 


0.742 


73M08C00 


0.508 


72N40CS4 


0.693 


72M40C84 


0.769 


72N72C84 


0.794 


T3N08C84 


0.952 


72N72C84 


0.744 


73N4CC28 


0.508 


72M72C84 


0.694 


72W72C34 


0.769 


72N40C84 


0.797 


72N40C84 


0.954 


73N08CS4 


0.753 


73M7*2C56 


0.508 


73M08C84 


0.720 


73M08C84 


0.797 


73N08C8^ 


0.809 


72N72C84 


0.954 


73N40CS4 


0./79 


74N72CS4 


0.503 


73W72C84 


0.742 


73N40C84 


0.831 


T3N72C34 


0.858 


73N40C84 


0.960 


73N72CS4 


0.781 


75NO8C00 


0.508 


73M40CS4 


0.744 


73M72C84 


0.834 


T3N40C34 


0.S64 


74N08CS4 


0.964 


74N08C84 


0.791 


75N40CS6 


0.508 


74M08CS4 


0.764 


74N08C84 


0.853 


74N08C84 


0.865 


"3H72C34 


0.968 


75N08CS4 


0.805 


T5N40CO0 


0.509 


74M40C34 


0.782 


74N40C84 


0.875 


75N08C34 


0.893 


74N40C34 


0.963 


7*N4CCS4 


0.807 


72M72C84 


0.510 


74M72C84 


0.788 


75M08C84 


0.875 


74^40C84 


0.905 


74N72CS4 


0.969 


74N72C84 


0.811 


74NC8C84 


0.510 


75N08CS4 


0.791 


74M72C84 


0.877 


74N72C84 


0.912 


75N02C84 


0.970 


75N72CS4 


0.820 


75N08C56 


0.510 


75N72CS4 


0.807 


75M72C84 


0.895 


75N72C84 


0.923 


75H40CS4 


0.973 


75N40CS4 


0.823 


72K72C00 


0.513 


75N'0CS4 


0.814 


75M40C84 


0.901 


75N40C34 


0.931 


75H72C34 


0.981 


* T: the 


numoer of 


treatments 


. H: the 


nuroer of s 


uojects 


oer treatment, C: the 


corre' .ion coeff 


c\ ent . 
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